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Abstract 

Purpose: Pterygium is a wing shaped fibrovascular growth on the ocular surface, characterized by fibrosis, angiogenesis, 
extracellular matrix remodeling, and inflammatory infiltrates. Epidemiologic studies have linked pterygium formation to 
various chronic inflammatory conditions, such as ultraviolet radiation, sawdust exposure, and dry eye disease. The purpose 
of this study is to identify proteins that are differentially expressed in primary pterygium by using a combination of gene 
microarray and proteomic platforms. 

Methods: Paired pterygium and uninvolved conjunctiva tissues of four patients were evaluated for differences in global 
gene transcript levels using a genechip microarray. Proteins extracted from another four pairs of tissues were quantified by 
iTRAQ approach. Western blot and immunofluorescent staining on additional patients were used to validate dysregulated 
protein expression obtained from microarray and proteomics data. In addition, primary conjunctival fibroblasts were treated 
with recombinant S100A8, S100A9 or both. Transcript level changes of a panel of potential target genes were evaluated by 
real time-PCR. 

Results:Tbe following were up-regulated at both protein and transcript levels S100 A8 and A9, aldehyde dehydrogenase 3 
family, memberl (ALDH3A1) and vimentin (VIM). Conversely, serpin peptidase inhibitor clade A member 1 (SERPINA1) and 
transferrin (TF) were down-regulated. Upon adding S100A8, S100A9 or both, the inflammatory chemokine CXCL1, matrix 
proteins vimentin, biglycan, and gelsolin, as well as annexin-A2, thymosin-pH, chymase (CMA1), member of Ras oncogene 
family RAB10 and SERPINA1 were found to be up-regulated. 

Conclusions: We identified 3 up-regulated and 2 down-regulated proteins by using a stringent approach comparing 
microarray and proteomic data. On stimulating cells with S100A8/9, a repertoire of key genes found to be up-regulated in 
pterygium tissue, were induced in these cells. S100A8/9 may be an upstream trigger for inflammation and other disease 
pathways in pterygium. 



Citation: Hou A, Lan W, Law KP, Khoo SCJ, Tin MQ, et al. (201 4) Evaluation of Global Differential Gene and Protein Expression in Primary Pterygium: SI 00A8 and 
S100A9 as Possible Drivers of a Signaling Network. PLoS ONE 9(5): e97402. doi:10.1371/journal.pone.0097402 

Editor: William B. Coleman, University of North Carolina School of Medicine, United States of America 

Received January 16, 2014; Accepted April 16, 2014; Published May 13, 2014 

Copyright: © 201 4 Hou et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This research is supported by the Singapore National Research Foundation under its NMRC/CSA/045/201 2, NMRC/BNIG/1 0075/201 2 and administered 
by the Singapore Ministry of Health's National Medical Research Council. The funders had no role in study design, data collection and analysis, decision to publish, 
or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: hou.aihua@seri.com. sg (AH); Louis.tong.h.t@snec.com.sg (LT) 
9 These authors contributed equally to this work. 



Introduction 

Pterygium is an ocular surface disease which characterized by 
abnormal epithelial and fibrovascular proliferation, invasion, and 
matrix remodeling. The lesion leads to astigmatism, ocular 
irritation, dryness, and reduced vision once it involves the visual 
axis. Despite its prevalence, the exact pathogenesis of pterygium 
remains unknown. Evidence from epidemiology and histopathol- 
ogy studies support the concept that pterygium arises as a result of 
ultraviolet (UV) radiation-damage [1-4]. UV radiation induces 
damage in DNA and other cellular substrates [5-7]. In addition, 



several studies have pointed to an immune response-inflammatory 
component in pterygium pathology [8-11]. 

Previously, microarray techniques have been used to compare 
the genome wide gene expression profiles between primary 
pterygium and normal conjunctiva tissues [12-14]. Extracellular- 
matrix-related, proinflammatory, angiogenic, fibrogenic, and 
oncogenic genes were found to be dysregulated in primary 
pterygium [13]. About 1724 genes were significantly different 
between primary and recurrent pterygium by genome-wide 
microarray scanning [15]. Our group reported 51 genes with at 
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least 2 fold changes between primary and recurrent pterygium 
[14]. More recently, proteomic techniques have been applied in 
pterygium research. Using two-dimensional electrophoresis and 
mass spectrometry (2D-MS), Kanamoto et al. identified 1 1 
differentially expressed proteins between primary and recurrent 
pterygium [16]. Using a similar assay technique, peroxiredoxin 2 
was found to be detected in primary pterygium tissue [1 7] . 

The S 1 00A8 and S 1 00A9 are calcium binding secreted proteins 
that have been previously identified to be up-regulated in various 
ocular surface conditions including pterygium, dry eye, meibo- 
mian gland dysfunction and corneal neovascularisation [18—21]. 
In living tissues, these two proteins may exist as homodimers or 
heterodimers and play a big role in regulation of oxidative damage 
[22-28] as well as initiation of immune responses [24,29-32]. 
Because of their possible role as a master regulator of immune 
response and inflammation, S100A8/S100A9 deserve further 
study. The relationship of S 1 00A8 and S 1 00A9 and downstream 
inflammatory molecules has not been evaluated in the context of 
the ocular surface. 

Since there are advantages to an unbiased molecular profiling 
method in a disease with uncertain etiology, we aimed to identify 
dysregulated molecules in primary pterygium tissue compared to 
normal uninvolved conjunctiva by proteomics evaluation and 
comparing findings with gene microarray analysis. We further 
investigated the ability of S100A8, S100A9 and S100A8/A9 to 
regulate the expression of a panel of genes involved in pterygium 
pathology using primary cell cultures. 

Methods and Materials 

1. Patients and Ethics Statement 

Primary nasal pterygium and uninvolved conjunctival tissues 
from the superior temporal quadrant of the same eye were 
collected and stored immediately in RNAfafer RNA Stabilization 
Reagent (Qiagen) for RNA extraction, or immediately frozen in 
liquid nitrogen for protein extraction. All procedures adhered to 
the tenets of the Declaration of Helsinki and were approved by the 
Institutional Review Board of the Singapore Eye Research 
Institute. Written informed consent was obtained from all 
participating patients. 

2. Microarray Gene Expression Evaluation 

RNA was extracted from 4 pairs of tissues using the following 
protocol: Tissues were minced into small pieces, and ground into 
powder using pestie and mortar with liquid nitrogen. The tissue 
powder was transferred into TRIzol (Invitrogen) and passed 
through 22G needles several times (Becton Dickenson). Chloro- 
form was added to each homogenized sample, mixed well and 
centrifuged at 13,000 rpm for 15 min at 4°C. The aqueous layer 
was transferred to a new tube and an equal volume of 70% ethanol 
added. Subsequendy the samples were purified using RNeasy 
RNA extraction kit (Qiagen) according to manufacturer's proto- 
col. RNA concentration and quality were assessed using the ND- 
1000 spectrophotometer (Thermo Scientific). 

The microarray experiment was performed on the Agilent 
Sureprint G3 human GE 8x60 K (ID: G685 1-605 10) single- 
colour chip by a commercial service (Molecular Genomics Pte Ltd, 
Singapore) according to the Agilent protocol. In summary, one- 
hundred nanograms of total RNA was labeled with Low Input 
Quick Amp Labeling Kit (One-Colour Microarray), then 
converted to double-stranded cDNA using oligo-dT primers and 
transcribed by T7 RNA polymerase to produce cyanine 3-CTP 
labeled cRNA. Sixty nanograms of labeled cRNA was hybridized 
onto Agilent SurePrint G3 human GE 8x60 K microarray for 



17 hrs at 65°C, 10 rpm in Agilent hybridization oven. After 
hybridization, the microarray slide was washed and scanned on 
Agilent High Resolution Microarray scanner. Raw signal data was 
extracted with Agilent Feature Extraction Software (V10.7.1.1), 
and analyzed by GeneSpring 12.5 (Agilent, CA). 

The threshold of raw signals was set to 1.0 and data was 
transformed to logbase 2 values. Per chip normalization, 
performed by baseline transformation to the median of all 
samples, was performed by percentile shift to 75 th percentile. A 
box plot of normalized data for the 8 samples showed a similar 
distribution. Subsequent quality control(QC) on the probeset to 
remove faulty or suspicious probes was done by filtering on flags 
(probes flagged as 'not detected' and 'compromised') and filtering 
on expression (removal of probes with very low raw signal intensity 
(<20)). This QC processing identified 27714 positive signals. The 
microarray data reported in this study had been deposited in 
NCBI Gene Expression Omnibus (GEO, http://www.ncbi.nlm. 
nih.gov/geo) with GEO series record number GSE51995. 

3. ITRAQ-LC-MS/MS Proteomics 

Protein extraction and quantification. Four pairs of tissues 
were weighed and washed with phosphate buffered saline (PBS), 
and then put into 100 \l\ of lysis buffer (75 mM NaCl, 1 mM 
EDTA, 0.2% CHAPS, 50 mM Sodium Fluoride, 1 mM sodium 
orthovanadate and 1 x protease inhibitors cocktail prepared in 
PBS). Ten freeze-thaw cycles involving snap-freezing with liquid 
nitrogen and thawing on ice water were subsequently performed. 
Samples were then ground using a mini grinder, and 400 JLll of the 
lysis buffer was added. Ground samples were further homogenized 
via sonication. Triton-X (0.2%) and IGEPAL (0.2%) were added 
after sonication. Homogenized samples were spun at 16,000 g for 
20 minutes at 4°C, the supernatant was concentrated to 
approximately 50 nL by a 0.5 ml 3 KDa ultra-centrifugal filter 
(Amicon, USA) and the protein concentration were measured by 
Bradford assay. 

Mass spectrometry and data analysis. Samples were 
reduced, alkylated, digested and labeled with iTRAQ following 
the manufacturer's instruction. The combined samples were 
cleaned and fractionated with strong cation exchange (SCX) 
chromatography on a PolySULFOETHYL A column (200x4.6- 
mm, 5 Urn, 200 A; PolyLC Inc., Columbia, MD, USA). Forty-one 
fractions were collected and they were subsequendy pooled to 15 
samples for mass spectrometry analysis. The collected fractions 
were first reduced volume and reconstituted in 1 % formic acid/ 
5% acetonitrile for desalting by solid phase extraction. C18 
purification was performed on Pierce C18 Spin Columns (Thermo 
Scientific). The desalted samples were lyophilized and reconsti- 
tuted in 20 ni of 0.1% formic acid/2% acetonitrile. LC-MS/MS 
analysis was performed on a QStar Elite quadruple time of flight 
mass spectrometer system (AB Sciex) coupled to a Tempo nano- 
LC (AB Sciex) with autosampler. 

Database search and relative analysis was performed on 
ProteinPilot™ software (AB Sciex) version 4.0.8085, rev 148085 
with the Paragon Algorithm™ version 4.0.0.0, rev 148083 using 
UniProtKB/Swiss-Prot database release 20 1 3_1 2 (containing 
49,243,530 entries) against the human genome (20274 sequences) 
[33-34]. Database search was performed by setting cysteine 
modification by MMTS as a fixed modification [35] and tryspin as 
the digesting enzyme. None of the special factors was selected. The 
unused Protscore (conf) in the software was set to 1.3 (95% 
confidence). Protein level changes were compared between paired 
samples of pterygium and conjunctiva, and iTRAQ fold change of 
>1.3 was considered to be up-regulated in pterygium, while a 
ratio of <0.77 was considered down-regulated. This is based on 
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fact that the systematic variation of our analytical system has been 
determined consistently to be within 30% across several studies 
[36-41]. 

4. Primary Pterygium and Conjunctival Fibroblast Cell 
Culture 

Primary pterygium and conjunctiva fibroblasts were cultured as 
previously reported [42-43]. Briefly, pterygium and conjunctiva 
tissues were immediately submerged in 1 ml of DMEM/F12 
(Gibco) medium containing 1 x antibiotic-antimycotic (1 x A/ A) 
(Invitrogen, USA) and transferred to lab on ice. Culture plates 
(60 mm) were coated with 1 ml of fetal bovine serum (FBS) 
(Gibco) for 10 min at 37°C. Each specimen was cut into pieces of 
approximately 1-2 mm and placed onto the FBS-coated plates. 
Sixty microliters of FBS with lx A/ A was dropped over each 
piece. Explants were incubated overnight at 37°C with 5% C0 2 . 
The next morning, 120 uj of DMEM/F12 with 10% FBS and 1 x 
A/A were dropped over each explant, and 400 jLLl more of the 
same medium was added later in the evening. About 3-5 days 
later, cells started to migrate from explants and the volume of 
medium was increased to 2 ml and changed every other day. 
Fibroblasts were then sub-cultured at day 18~21 using TrypL 
Express (Life Technologies). Cells from passage 1 to passage 3 
were used for experiments. Cell population homogeneity was 
observed as previously described [44]. 

5. Western Blot 

Western blots were performed as described [44]. In summary, 
total protein (30 (lg) from primary cultured pterygium and 
conjunctiva fibroblast cells were electrophoresed in 10—12% 
SDS-polyacrylamide gels, and then transferred to polyvinylidene 
difluoride (PVDF) membranes at 4°C overnight. The membranes 
were subsequently incubated with primary antibodies to 
ALDH3A1 (Abeam, UK, Cat. No. AB76976, 1:1000) and 
GAPDH (Abeam, UK, AB37168, 1:1000) for 2 hrs at room 
temperature. Membranes were then incubated with secondary 
antibodies conjugated to horseradish peroxidase for 1 hr at room 
temperature. Immunoreactivity was visualized with chemilumi- 
nescence substrate (SuperSignal West Pico or West Dura; Pierce). 

6. Immunofluorescent Staining 

Pterygium and conjunctival tissues were sectioned with a 
Microm HM550 cryostat (Microm, Walldorf, Germany) at 8 |Xm 
thickness. Sections were fixed with 4% paraformaldehyde for 
20 min, washed with PBS, permeated in PBS containing 0.15% 
Triton X-100 for 15 min, blocked with 4% BSA in PBS containing 
0.1% Triton X-100 (Sigma) for 1 hr, and then incubated with 
ALDH3A1 antibody (Abeam, UK, Cat. No. AB76976, 1:200 
dilution in blocking solution) at 4°C overnight. After washing with 
PBS containing 0. 1 % Tween-20, the sections were incubated with 
Alexa Fluor 488-conjugated secondary antibody at room temper- 
ature for 45 min. Slides were then mounted with UltraCruz 
Mounting Medium (Santa Cruz). For negative controls, non- 
immune reactive serum was used in place of the specific primary 
antibody. Sections were observed and imaged using a Zeiss 
Axioplan 2 fluorescence microscope (Zeiss, Oberkochen, Ger- 
many). 

7. S100A8, S100A9 and S100A8/A9 Treatment of 
Conjunctiva Fibroblast Cells 

The maximum non-toxic dose of recombinant S100A8, S100A9 
and combined S100A8/S100A9 for primary conjunctiva fibro- 
blasts treatment were pre-determined to be 40, 25 and 20 |Ig/ ml 



respectively. Primary conjunctiva fibroblast cells were seeded 
(200,000 per well) in a 12-well culture plate (Greiner bio-one) and 
cultured in serum-free DMEM/F12 medium (Gibco) for 24 hrs. 
Cells were washed with PBS and treated with S100A8, S100A9 
and S100A8/A9 proteins in fresh media for 6 or 24 hrs. Control 
cells were given fresh medium without recombinant protein. 

8. Real Time PCR 

RNA was extracted with RNeasy mini kit (Qiagen) according to 
the manufacturer's instructions. Reverse transcription was per- 
formed using Superscript III first strand synthesis supermix 
(Invitrogen) according to manufacturer's protocol. Real time- 
relative quantitative polymerase chain reaction (qPCR) was 
performed with SYBR Green qPCR Mastermix (Roche). Primers 
used for real time PCR are listed in Table 1. GAPDH was used as 
the endogenous control gene. PCR cycles were performed on the 
Roche LightCycler 480 with the following cycle conditions: 
10 min at 95°C, 45 cycles of 10 s at 95°C, 15 s at 60°C and 
20 s at 72°C, followed by melting curve analysis. The delta-delta 
Ct method was used for data analysis and experimental 
reproducibility was confirmed using 3 biological replicates. 
Statistical significance was evaluated at p<0.05. Pooled standard 
deviations of the fold changes were calculated for 3 independent 
experiments which used cells cultured from 3 different patients. 

Results 

1. A Panel of Genes Found to be Dysregulated at both 
Transcript and Protein Level 

Fold change analysis was performed using the microarray data 
from the 4 paired pterygium and conjunctiva tissues, a list of 
differentially expressed genes (fold change ^2, p<0.05) was 
delineated (GEO series record number GSE51995). Some of the 
significandy dysregulated transcripts in pterygium are shown in 
Table 2. Other up-regulated transcripts (Table SI) include 
chemokine ligand CXCL1, biglycan (BGN), and gelsolin (GSN), 
annexin-A2 (ANXA2), thymosin-fJ4 (TMSB4X) and member of 
Ras oncogene family (RAB10), whereas the transcripts down- 
regulated in pterygium were those of Serpin peptidase inhibitor, 
clade A member l(SERPINAl) and chymase (CMA1). At least 3 
out of 4 pairs of the transcripts were regulated in the direction 
shown although p<0.05 was only for GSN and SERPINA1. We 
confirmed that in the dataset (http:/ /www.ncbi.nlm.nih.gov/geo/, 
series accession GSE2513) from our previous study [14] these 
transcripts were also dysregulated in the direction shown (Table 
SI). In the previous dataset, the p values were <0.05 for 3 of these 
genes: BGN, GSN and TMSB4X, and these were 0.003, 0.033 
and 0.014 respectively. In the current dataset, however, none of 
the evaluated genes, passed the Benjamini-Hochberg False 
discovery rate correction (p<0.05) when the paired student's T- 
test was used, suggesting heterogeneity amongst different biolog- 
ical samples. 

The iTRAO_ LC-MS/MS experimental design is shown in 
Figure SI. The complete raw data of the iTRAQ, analyses are 
shown in Table S2 and S3 respectively. Only proteins identified 
with 95% confidence were included. After relative quantification 
analysis, all the proteins found to be significandy different between 
pterygium and conjunctiva at least in 3 pairs of tissues are shown 
in Table S4. Inter-patient heterogeneity of protein expression was 
also observed across many proteins. When the list of dysregulated 
proteins in proteomics was compared against that from the 
microarray data, a panel of overlapping genes was generated and 
summarized in Table 2. Although not all functionally important 



PLOS ONE | www.plosone.org 



3 



May 2014 | Volume 9 | Issue 5 | e97402 



S100A8 and S100A9: Possible Drivers of a Signaling Network in Pterygium 



Table 1. Primers used in Real Time PCR. 







Gene 
Symbol 


Gene Name 


NCBI accession number 


Primer Sequence 


TMSB4X 


Thymosin beta 4, X-linked 


NM_021 109.3 


F: AGACCAGACTTCGCTCGTA 
K: CTGCTTGCTTCTCCTGTT 


CXCL1 


Chemokine(C-X-C motif) ligand 1 


NM_001511.3 


F: CCAGACCCGCCTGCT 
R:CCTCCTCCCTTCTGGTCAGTT 


ANXA2 


Annexin A2 


NM_001 002857.1 


F: CTCTACACCCCCAAGTGCAT 
R: TCAGTGCTGATGCAAGTTCC 


SERPINA1 


Serpin peptidase inhibitor, clade A member 1 


NM_000295.4 


F: GTCAAGGACACCGAGGAAGA 
R: TATTTCATCAGCAGCACCCA 


CMA1 


Chymase 1 


NM_001 836.3 


F: CTGGGGAAGAACAGGTGTGT 
R: 1 1 1 1 GTCTTCCTGGG ATTGC 


GSN 


Gelsolin 


NM_000177.4 


F: TGGTGGTGCAGAGACTCTTCC 
R: CTTTCATACCGATTGCTGTT 


BGN 


Biglycan 


NM_001711.4 


F: CCTTCGACCAGTCCTCCCTTC 
R: GGGACAGGCGAAGCCAGGTTC 


VIM 


Vimentin 


NM_003380.3 


F: CCAAACTTTTCCTCCCTGAACC 
R:GTGATGCTGAGAAGTTTCGTTGA 


RAB10 


Ras-related protein Rab-10 


NM_016131.4 


F: CACTTCTTCCTTGCGTCTCC 
R: GAAGGCATCTGGGACACATT 


GADPH 


Glyceraldehyde 3-phosphate dehydrogenase 


NC_000012.11 


F: CCCCACACACATGCACTTACC 
R: CCTACTCCCAGGGCTTTGATT 



doi:1 0.1 371 /journal.pone.0097402.t001 



molecules are dysregulated at both the transcript and protein 
levels, we decided to report and emphasize these for simplicity. 

Aldehyde dehydrogenase 3 family, memberl (ALDH3A1), 
vimentin (VIM), S100A8 and S100A9 were found to be up- 
regulated in primary pterygium compared to healthy conjunctiva 
control. SERPINA1 and transferrin (TF) were down-regulated at 
both protein and transcript levels in pterygium (Table 2). 

2. Validation of Over-expression of ALDH3A1 by Western 
Blot and Immunofluorescent Staining 

On immunofluorescent staining in pterygium and conjunctiva 
tissues (Figure 1A, left and middle columns), ALDH3A1 was 
found to be present in both conjunctival and pterygium 
epithelium, with a stronger expression in pterygium than in 
conjunctiva. The localization of the ALDH3A1 was in the 
cytoplasm and cell periphery of the epithelium, and present in 
all layers of the epithelium. Some fibroblast staining was visualized 



in the stroma of both pterygium and conjunctiva but comparison 
of staining intensity was difficult (data not shown). No fluorescent 
signals were detected in the negative controls (Figure 1A, right 
column). Levels of ALDH3A1 in primary cultured pterygium and 
conjunctiva fibroblast cells and epithelial cells were analyzed by 
western blot. Results showed that ALDH3A1 was more highly 
expressed in pterygium fibroblast cells than conjunctiva fibroblast 
cells (Figure IB). However, we were not able to show similar up- 
regulation in cultured epithelial cells (data not shown). 

3. S100A8, S100A9 and S100A8/S100A9 Up-regulate 
Disease-inducing Genes 

Previously, our group has reported that S100A8 and S100A9 
were localized in the superficial layer of both pterygium and 
normal conjunctiva epithelium, with higher levels in pterygium 
than uninvolved conjunctiva [21]. To test our hypothesis that 
S100A8/S100A9 are master regulator of many other dysregulated 



Table 2. Genes dysregulated in both transcriptomics and proteomics. 



Microarray Proteomics 



Average iTRAQ fold change* 




Gene Symbol 


Description 


P 


Fold change 


(p<0.05) 




ALDH3A1 


aldehyde dehydrogenase 3 family, memberl 


0.01 


2.9 


8.24 


Up regulated 


VIM 


vimentin 


0.45 


1.1 


4.64 




S100A8 


S100 calcium binding protein A8 


0.47 


1.2 


3.98 




S100A9 


S100 calcium binding protein A9 


0.64 


1.2 


5.93 


Down regulated 


SERPINA1 


serpin peptidase inhibitor, clade A member 1 


0.01 


3 


0.29 




TF 


transferrin 


0.24 


1.5 


0.48 



*averaging obtained in duplicated analyses of 4 pairs of tissues. 
doi:1 0.1 371 /joumal.pone.0097402.t002 
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A 





ALDH3A1 
GAPDH 



Figure 1. ALDH3A1 expression in pterygium. A. Protein expression of ALDH3A1 in both conjunctiva and pterygium fibroblast cells was 
determined by immunoblotting with antibodies specific for ALDH3A1 and GAPDH. B. Immunofluorescent staining image showing presence of 
ALDH3A1 (green) in human conjunctiva and pterygium epithelium. Nuclei were stained with DAPI (blue) present in the mounting medium. Negative 
control (without primary antibody) was shown. All images were taken at 200 x magnification. C: conjunctiva; P: pterygium. 
doi:1 0.1 371 /journal.pone.0097402.g001 



genes, recombinant human S100A8, S100A9 or both were used to 
treat conjunctiva fibroblast cells. A panel of dysregulated genes in 
pterygium (mentioned above) was selected as possible downstream 
effectors of the S100 proteins. The transcript level changes of these 
genes in cultured conjunctival fibroblasts, induced by addition of 
recombinant S100A8, S100A9 or both were evaluated by qPCR. 

Inflammation-related genes TMSB4X, CXCL1 and ANXA2, 
were significantly increased upon addition of S100A8, S100A9 
singly or in combination (Figure 2). TMSB4X transcript was 
significantly up-regulated by any of these three interventions at 
6 hrs and continued to increase 24 hrs after treatment, while a 
significant up-regulation of TMSB4X by S100A9 was observed 
only at 24 hrs after treatment (Figure 2A-C). CXCL1 showed a 
marked fold change (of more than 700 folds) upon S100A8 and 
S100A9 treatment at 6 hrs, but this was less marked at 24 hrs. 
CXCL1 expression was not responsive to S100A8/A9 combina- 
tion at 6 hrs but showed an increase at 24 hrs (Figure 2D-F). 
ANXA2 transcript increased significantly at 6 hrs after any of 
these interventions, and the increased expression persisted through 
24 hrs, but this transcript was relatively more sensitive to the 
addition of the S100A8/A9 combination than single proteins. 

Extracellular matrix genes GSN, BGN and VIM increased 
upon S100A8, S100A9 or combined treatment in conjunctiva 
fibroblast cells (Figure 3), but the magnitude of change (1.5 to 20 



fold) was generally not as large as the transcripts for inflammation- 
related genes. 

The expression levels of enzyme-coding genes SERPINA1 and 
CMA1 were up-regulated by S100A9 (Figure 4B and 4E 
respectively) and combination of S 1 00 A8 and S 1 00A9 (Figure 4C 
and 4F), but not by S100A8 (Figure 4A and 4D). Intracellular 
trafficking gene RAB10 was up-regulated by S100A8 and S100A9 
at 6 hrs (Figure 4G and 4H), while the increase persisted 6 hrs 
through 24 hrs for the combination treatment (Figure 41). 

Discussion 

In this study we reported a list of differentially expressed genes 
by genome-wide microarray analysis using pairs of primary 
pterygium and autologous conjunctiva tissues, and identified 10 
dysregulated proteins by iTRAQ;LG MS/MS proteomics. Five 
genes were found to be dysregulated in pterygium tissue in the 
same directions using both approaches. In addition, we found that 
S100A8, S100A9 and combination could up-regulate genes 
related to inflammation, extracellular matrix and matrix enzymes 
in cultured primary human conjunctiva fibroblast cells. 

A microarray study conducted by John-Aryankalayi et al 
compared two primary and one recurrent pterygium with 
autologous conjunctiva tissues [13]. Previously, we performed a 
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Time after S100A8 treatment (hr) Time after S100A9 treatment (hr) Time after S100A8/A9 treatment (hr) 




Time after S100A8 treatment (hr) Time after S100A9 treatment (hr) Time after S100A8/A9 treatment (hr) 



Figure 2. Real-time PCR quantification of fold changes in TMSB4X, CXCL1 and ANXA2 transcripts in human primary conjunctiva 
fibroblasts after addition of S100A8 (left column, white bars), S100A9 (middle column, grey bars) and both (right column, black 
bars) recombination proteins. Cells were incubated with these proteins for 6 hrs and 24 hrs. The level of the transcript was normalized to the 
control cells. Height of the bars represents the mean of the relative fold changes. Error bars represent the standard deviation among 3 patient 
samples. Indicates p<0.05, ^indicates p<0.01 and ***p<0.001. 
doi:1 0.1 371 /journal.pone.0097402.g002 



microarray study with 4 pooled samples of conjunctiva tissues and 
10 un-pooled samples of primary pterygium tissues [14]. Different 
with all previous studies, our current microarray utilized 4 pairs of 
primary pterygium and conjunctiva tissues from the same eyes. 
The differentially regulated genes that were further evaluated in 
the current study were also found to be dysregulated in the same 
direction in the previous dataset uploaded to GEO (http:/ /www. 
ncbi.nlm.nih.gov/geo/, series accession GSE2513) [14]. 

To date, there have been only two published studies which 
identified dysregulated proteins in pterygium using proteomics 
[16-17]. Both published method utilized 2-D PAGE for protein 
separation and detection. While 2-D PAGE has powerful 
resolution, it is generally incapable of resolving proteins of extreme 
PI and molecular weight. Furthermore, sensitivity of the analytical 



system is largely dependent on the ability to visual protein spots on 
2-D gels. In contrast, this study used a more modern technology 
involving isobaric tagging for relative and absolute quantification 
(iTRAQ). Coupled to LC/MS/MS, this method is capable of 
analyzing up to 8 samples in a single analysis thus providing a 
more reliable comparative analysis by eliminating the reproduc- 
ibility issue encountered in 2-D PAGE. One of the 2-D PAGE 
based studies identified 1 1 proteins differentially expressed 
between recurrent and primary pterygium [16], but did not 
examine differences between primary pterygium and normal 
control tissue. The other study reported the over-expression of 
only 1 protein, peroxiredoxin 2, in primary pterygium [17]. 

The dysregulated genes identified were known to have 
biological functions relevant to inflammation or cell stress. 
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Figure 3. Real-time PCR quantification of GSN, BGN and VIM transcripts fold change in human conjunctiva fibroblast after S100A 
proteins treatment. Refer to Figure 2 for details. 
doi:1 0.1 371 /journal.pone.0097402.g003 



ALDH3A1 catalyzes the oxidation of a wide variety of endogenous 
and exogenous aldehydes to the corresponding carboxylic acid, 
and is associated with cell proliferation and drug metabolism [45] . 
This enzyme is a corneal crystallin and plays a critical role in 
protecting the cornea from UV-induced oxidative stress by direcdy 
absorbing UV-radiation [46-48], scavenging free radicals or 
producing NADPH [45]. A previous study had found ALDH3A1 
to be among the 25 most abundant transcripts in an un- 
normalized un-amplified pterygium cDNA library [13]. Over- 
expression of ALDH3A1 has also previously been found in tumor 
cells and increased ALDH3A1 expression can lead to cell 
proliferation [45]. The dysregulation of ALDH3A1 in pterygium 



tissue may be attributed to abnormal control of oxidative stress or 
proliferation. ALDH3A1 in cultured pterygium and conjunctiva 
epithelial cells showed no significant difference, most likely is 
related to differences in ALDH3A1 in cultured and native 
epithelium, the latter having matrix and fibroblast interactions. 

Vimentin is a type III intermediate filament, a major 
cytoskeletal component of mesenchymal cells, and present in cells 
at the focal areas of a surgical pterygium tissue [49] . It has been 
reported to be involved in corneal wound healing and fibrosis [50], 
which are processes that may be implicated in pterygium 
pathology. In fact, abnormal vimentin expressing cells were 
reported in pterygium tissue [51]. 
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Figure 4. Real-time PCR quantification of SERPINA1, CMA1 and RAB10 transcripts fold change in human conjunctiva fibroblast 
after S100A proteins treatment. Please refer to Figure 2 for details. 
doi:1 0.1 371 /journal.pone.0097402.g004 



The S100A proteins are a family of small calcium-binding 
dimeric proteins with pleiotropic functions, and the S100A8/A9 
heterodimer play an important role in inflammatory regulation 
[24,52]. S100A proteins are abundantly expressed on the ocular 
surface [18] and have been suggested to be involved in the 
pathogenesis of various ocular surface diseases, such as dry eye 
[20], corneal angiogenesis [19] and pterygium [20-21]. 

SERPINA1, also known as a 1 -proteinase inhibitor or 0(1- 
antitrypsin, is a protease inhibitor. Decreased levels of SERPI- 
NAlat both transcript and protein levels were noticed in 
circulating hematopoietic cells, induced by granulocyte colony- 
stimulating factor or chemotherapy [53]. TF is an iron-binding 



blood plasma glycoprotein that regulates the level of free iron in 
biological fluids [54]. 

Pterygium pathogenesis is associated with increased UV 
exposure [1-4] and S100A8 protein is known to be induced by 
UV radiation [55-56]. Functions of the genes induced by addition 
of S 100 A proteins may suggest that they propagate the pterygium 
disease. It is possible that the S100A proteins may serve as a 
master regulator to activate this network of genes (Figure 5). In 
our proposed hypothesis, SERPINA1 and ALDH3A1 may also be 
dysregulated under the influence of UV radiation, but may not be 
regulated through S100A8/A9. Nevertheless, Figure 5 shows 
how these pathways may be linked. For example, the ANXA2 is 
involved in resistance to UV-induced cell death [57]. Although we 
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Figure 5. Schematic showing how UV may trigger upstream molecular disturbances such as S100A8 and S100A9, ALDH3A1, and 
SERPINA2 which result in a cascade of events leading to pterygium progression. [ represents up-regulation in pterygium compared to 
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used fibroblasts in the study, other cell types such as dendritic cells, 
may have important responses to SI 00 A proteins in pterygium 
and should be studied further. 

A strength of the current study is that comparisons were done 
using un-pooled samples from individuals, avoiding noise related 
to inter-individual differences. As the amount of RNA extracted 
from small surgical samples of pterygium and conjunctiva was 
limited, earlier genome-wide gene expression studies required 
pooling of samples from different patients [12,14-15]. Another 
shortcoming of sample pooling is that some of the differentially 
expressed genes may be obscured. 

A major limitation of the protein studies is that we did not do a 
more definitive evaluation of small molecules like cytokines using a 
technique such as bead-based indirect immunofluorescent assay, 
as such we may not be able to verify the expression of some 
chemokines and cytokines. Another limitation of this study is that 
pterygium and conjunctiva tissues have different epithelium/ 
stroma ratio which definitely introduce variations. This should be 
considered when interpreting the list of dysregulated molecules. 

In conclusion, several genes were dysregulated at the transcript 
and protein levels in pterygium. The dysregulated molecules were 
consistent with previous experiments reported. The S100A8 and 
S100A9 have been found to be up-regulated in pterygium tissue. 
These may control the progression of pterygium via a network of 
other molecules which have also been found to be dysregulated in 
pterygium. 
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Figure SI Four pairs of pterygium and conjunctiva samples 
were used for two 8-plex iTRAQ. LC-MS/MS. El and El 
represented two 8-plex experiments. P: pterygium; C: conjunctiva, 
i and ii: each pair of samples was duplicated. 
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